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Abstract 

We discuss the evolution of the mixed phase at RHIG and SPS within boostinvari- 
ant hydrodynamics. In addition to the hydrodynamical expansion, we also consider 
evaporation of particles off the surface of the fluid. The back-reaction of this evapo¬ 
ration process on the dynamics of the fluid shortens the lifetime of the mixed phase. 
In our model this lifetime of the mixed phase is < 12 /m/c in Au -|- Au at RHIG and 
< 6.5 frujc in Pb -|- Pb at SPS, even in the limit of vanishing transverse expansion 
velocity. Strong separation of strangeness occurs, especially in events (or at rapidities) 
with relatively high initial net baryon and strangeness number, enhancing the multi¬ 
plicity of MEMOs (multiply strange nuclear clusters). If antiquarks and antibaryons 
reach saturation in the course of the pure QGP or mixed phase, we find that at RHIG 
the ratio of antideuterons to deuterons may exceed 0.3 and even '^He/^He > 0.1. In 
S + Au at SPS we find only N/N ~ 0.1. Due to fluctuations, at RHIG even negative 
baryon number at midrapidity is possible in individual events, so that the antibaryon 
and antibaryon-cluster yields exceed those of the corresponding baryons and clusters. 


Work supported bp BMBF, DFG, GSL 
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1 Introduction 


The initial (pre-equilibrium) stage of heavy-ion collisions at BNL-RHIC energies has been 
intensively studied during the last years within various models, see e.g. refs. [||, |^, |^. As 
a consequence of the huge number of (sea) partons in the inital state and the large cross 
sections (especially for gg) already after a very short (proper) time r* ~ 0.1 — 0.2 fm/c the 
parton momentum space distributions at midrapidity are isotropic. Their widths correspond 
to a very high temperature on the order of T, ~ 500 MeV. On the other hand, the approach 
to chemical equilibrium probably takes somewhat longer. Nevertheless, it appears reason¬ 
able to assume that at midrapidity a quark-gluon plasma (QGP) in thermal and chemical 
equilibrium emerges after a few fm/c. This QGP subsequently expands and cools. If equi¬ 
librium is maintained during the expansion and if the phase transition to the hnally observed 
hadrons is of hrst order (with a sufficiently large latent heat), the system will spend a large 
part of its lifetime in the mixed phase BBB- 

In this paper we discuss a scenario for the hydrodynamical evolution of the mixed phase, 
also including particle evaporation off the surface of the system and the back-reaction of this 
time-dependent freeze-out. We also account for net strangeness conservation to compute 
the multiplicities of strange hadrons and -clusters. General aspects and thermodynamic 
properties of strange matter are discussed in ^ . 

We do not consider the evolution of the pure QGP phase since most calculations of the 
initial conditions indicate that at early times r < 1 — 3 fm/c the parton gas might not be in 
chemical equilibrium, even though it is thermalized Also, the lifetime of the mixed 

phase at RHIG exceeds that of the pure QGP by roughly a factor of 3. Thus, probably many 
more hadrons are evaporated from the mixed phase. Anyhow, the emission of hadrons from 
the pure QGP phase m would further increase the evaporation rates (cf. section |^) and the 
effects discussed here would be even stronger. 

We consider the following dynamical picture. At (proper) time tm a chemically and 
thermally equilibrated system is formed, which is assumed to be in the mixed phase. The 
mixed phase consists of a quark-gluon plasma modelled by a MIT bagmodel equation of 
state (mass of s-quark: = 150 MeV), and a hadron gas that includes all well-established 

(strange and non-strange) hadrons up to masses of 2 GeV (the list is given in [|^]). For the 
hadron gas Hagedorn’s eigenvolume correction is applied [Q. We employ a bag parameter 
of B = 380 MeV/fm^, leading to a mixed phase temperature (at small fig and /r^) of 
Tc ~ 160 MeV. The same equation of state was used in to study hadron ratios, cluster 
abundancies, and strangeness distillation in S' -|- Au reactions at GERN-SPS within a similar 
model B |T^, see also p. However, in contrast to the model presented here, spherical 


symmetry was assumed and boostinvariant hydrodynamical expansion was not taken into 
account. Also, the evaporation rates (cf. section Q) are different in the two models. 

Initially the fraction of volume occupied by quarks and gluons equals one, \{tm) = 
yQGP jytot _ Hydrodynamic expansion and evaporation of hadrons from the surface set 
in, leading to a decrease of the total energy-, entropy-, and net baryon density. Except 
for the permanent decoupling of hadrons at the boundary to vacuum, the mixed phase is 
assumed to remain in equilibrium while hadronization proceeds. 
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At some (proper) time th the quark-gluon content of the mixed phase (i.e., A) vanishes 
and all energy etc. has been converted into hadrons. In principle, the hydrodynamic evolution 
continues until the decoupling stage is reached. However, experimental data for heavy-ion 
collisions at AGS and SPS indicate that the chemical potentials and the temperature at the 
time of decoupling are close to the phase boundary |^], i.e. that the hydrodynamic expansion 
of the hadronic phase can be neglected. We therefore decouple all remaining hadrons from the 
hydrodynamic flow immediately after the conversion of mixed phase matter into hadrons is 
complete. Note however, that especially the multiplicity of antibaryons and anticlusters may 
depend strongly on the volume and the lifetime of the hadronic phase (due to annihilation) 
jl^, [l4| , even if it is very short. 

We also employ the same dynamical picture (with different initial conditions, of course) 
for S + Au and Ph + Ph collisions at SPS, although the assumptions of boostinvariant expan¬ 
sion and full strangeness saturation might only be moderately fulhlled at this energy. Nev¬ 
ertheless, by comparing our results to those of ref. (evaporation without boostinvariant 
expansion), we can give upper and lower limits for the time-scale of chemical freeze-out. 


2 Hydrodynamic expansion including particle evapo¬ 
ration 


The dynamical evolution of the mixed phase is described within a simple hydrodynamical 
model which is consistent with the inside-outside cascade picture at high energies [|^. In 
this model, the thermodynamic state of the system is dehned along proper time hyperbolas 
r = \/P — z^. The flow-rapidity r] is assumed to be equal to the space-time rapidity, i.e. 
the (longitudinal) flow velocity at the space-time point {t,z) equals ny = zjt. The space- 
time volume element is d^x = d^XTTdrdr]. In the case of a perfect fluid (i.e. disregarding 
viscosity and thermal conductivity), the rapidity density of each quantity corresponding to 
a conserved current is (proper) time independent. 


d dAjr) 

dr dr] 


( 1 ) 


In the case without explicit particle evaporation, this equation holds for the (total) entropy 
S, the net baryon number Nb, and the net strangeness Ng. The energy per rapidity interval 
decreases due to the work performed by the expansion according to 


d dE{T) 
dr dr] 


= -p 


d dV{T) 
dr dr] 


= —tiRbP 


( 2 ) 


where p denotes the pressure. For simplicity, we consider only longitudinal expansion, i.e. 
dPr/dr = 0, assuming that collective transverse expansion of the mixed phase is slow. The 
transverse shock wave, which can additionally produce up to ~ 10% of entropy (for almost 
net baryon free matter) p, |^, therefore is absent. A microscopic calculation of entropy 
production in the phase transition to hadronic matter can be found in ref. fl^. In our case. 
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this additional entropy can however be thonght of being absorbed in the initial entropy which 
has comparable nncertainties (see next section). 

In addition to the hydrodynamic expansion we consider particle evaporation off the snr- 
face of the midrapidity cylinder with transverse radius Rt- For the evaporation rate in the 
local restframe of the fluid we estimate 


drdr] 


O 


H 


f (■p\PT-rT ^ . 

— fi{E) -0 [pt ■ Vt) 

E Tt 




O^n 


vr 



(3) 


where = 2ttRtt denotes the surface of the expanding cylinder, is the density of the 
hadron species i under consideration. 

The evaporated particles carry away entropy, baryon number, strangeness, and energy: 


dS^^ _ 

drdr] ^ * drdr] ’ 

drdr] ^ drdr] ’ 

diVf _ _ dNf^ 

drdp ^ drdr] 


(4) 

(5) 

( 6 ) 


Si = sf /nf, Ms,! = nfi/nf, and AfB,i = denote the entropy, strangeness, and 

baryon number per particle of the hadron species i (for the given values of ps, I^-b, T), 
respectively. For the pions and temperatures around 160 MeV we compute ~ 0.75. 

For all other particles, we employ 



which is close to the non-relativistic limit 


{vt)7 = 



(7) 

( 8 ) 


By numerical integration of eq. (^ we have found that in the relevant temperature range (|^) 
gives a better approximation of the exact value than (H), especially for kaons and nucleons. 
The evaporated energy could be computed employing an expression similar to eqs. (@)-(||), 


dE^'^ 

drdrj 




dNf 

drdp 


(9) 


with Si = ef /nf denoting the energy per hadron of species i. Eq. (||) thus assumes that the 
energy per particle in the free-streaming regime is equal to that in the thermal ensemble. 
However, this would mean that the evaporation of hadrons performs mechanical work, i.e.. 


- pdV^^ = dE^^ - TdS^^ - pidN^^ > 0 


( 10 ) 
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The volume of the hadron gas simply shrinks until it disappears, while the volume of the QGP 
stays constant. Thus, eqs. 0-© do not lead to conversion of QGP into hadronic matter, 
unless transverse expansion is also taken into account (by adding a term —2'KRTTpdRT/dT 
on the right-hand-side of eq. (^)). We have found that for moderate transverse expansion 
velocities, dRx/dT 0.2, the shrinking of the transverse radius due to evaporation and the 
increase due to expansion cancel. The resulting evolution and lifetime of the mixed phase 
are very similar to that described below. 

We therefore assume that evaporation causes only heat transfer (TdS'®’' < 0) and chemical 
work {pdN^'" < 0) but does not perform mechanical work. In this case the evaporated energy 
is determined by the hrst law of thermodynamics (with = 0), 


dE^'^ _ dS^^ dNl^ dN^^ 

drdrj drdr] drdr] dr dr] 

The time evolution of the thermodynamic quantities due 
hydrodynamic expansion is 


( 11 ) 


to particle evaporation and 


dA 

drdrj 

dE 

drdrj 


dA^^ 

drdrj 


(A = iV„ Nb) 


dE^^ 

drdrj 


— ttR^p 


( 12 ) 

(13) 


These equations generalize (Q,^ to account for the back-reaction of the freeze-out on the 
dynamics of the fluid. Note that there are two distinct time scales in this model, one given 
by the (longitudinal) hydrodynamic expansion and the other by the evaporation process. 

The preequilibrium parton evolution yields the initial values of the entropy, the net baryon 
number, and the net strangeness, i.e. dSijM)/dp, dNB{rM)/dp, and dNgijM)/dp (see below). 
By virtue of the equation of state and the additional requirement that the system is in the 
mixed phase with a QGP-fraction of 100% (i.e., \{rM) = 1), the initial energy at midrapidity 
is also determined. 

yield dE/dp, dS/dp, dNB/dp, and dNg/dp. From 

we can compute 

jiq, jis, T, and the volume fraction of QGP A = /{V^ -|- 

The evaporated particles are assumed to stream freely to the detector. The measured 
hadron multiplicities are thus given by the evaporation rates (^, integrated over (proper) 
time. In addition, at the time rn where the hadronization is complete {X{rH) = 0), the entire 
remaining hadronic system is assumed to decouple instantaneously (i.e., on the hyperbola 
r = rn) from thermal and chemical equilibrium. 


At later times r > rM, eqs. (|T2 


these four quantities and Gibbs’ condition of phase equilibrium, = p^ 


3 Initial conditions 

On average over many events, the natural choice for the initial net strangeness (at midrapid¬ 
ity) is dNs/dp = 0. In contrast to heavy-ion reactions at fixed-target energies, the net baryon 
number at midrapidity is expected to be relatively small at RHIG. This is due to the large 
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number of secondary partons as compared to the number of valence quarks. Nevertheless, 
it is important to account for nonvanishing Nb, especially when studying the enrichment of 
the mixed phase with net strangeness. The models RQMD 1.07 FRITIOF 7.02 
and the PCM predict dNs/dr] = 20 — 35 in Au + Au at RHIC. 

The entropy per baryon created in the early evolution of the parton gas is S/Nb = 
200 in the PCM [^. In ref. S/Nb = 110 was estimated for the “final” (i.e. at the 
time where isentropic expansion sets in) entropy per net baryon by assuming perturbative 
QCD interactions (lowest order) between the partons in the initial state. However, in this 
calculation a transverse momentum cut-off (which is necessary to render the lowest-order 
pQCD cross sections hnite) of po = 2 GeV was used. As discussed in |^, this might be 
appropriate for LHC but for lower energies like RHIC the saturation of the produced gluons 
should be reached already for a smaller cut-off. Changing po from 2 GeV to 1 GeV increases 
the minijet cross sections and thus S/Nb to about 150 
assumed to be constant from SPS to RHIC, {S/Nb)^^^^ 


Also, the soft component was 
50. Multiplying this number by 
an ad-hoc factor of 2 increases the total S/Nb to about 150 for po = 2 GeV and to 220 for 
Po = 1 GeV |2|. 


The RQMD 1.07 model, based on color-string formation and decay, predicts a central 

|0 


pion multiplicity of dN.„/dr] = 1000, leading to St,/Nb = S.QNt^/Nb = 170 


If we 


assume that other particles (e.g. kaons) contribute about 10% — 20% to the total entropy, 
we hud a similar number for S/Nb as in the PCM. The microscopic phase-space model 
UrQMD 0, which also approaches RHIC from below assuming that color-string formation 
and subsequent fragmentation is the dominant mechanism for secondary production (as at 
SPS), predicts about 30 net baryons, 900 pions, and 80 —100 kaons at midrapidity in central 
(6 = 2 fm) Au + Au collisions at RHIC Thus, a value of S/Nb = 200 (at midrapidity) 
for an average central Au -\- Au event at RHIC seems reasonable. 

Considering tm (the initial time for the mixed phase), we first note that once the midra¬ 
pidity numbers for Nb and S are hxed, the corresponding coordinate space densities ub and 
s behave like 1/r. Thus, tm must not be chosen too small for then the system cannot be 
in the mixed phase initially (the equation pP*^^ = p^ cannot be fulhlled). Also, the initial 
transverse radius Rt, as given by the thermodynamic relations, depends on tm (through the 
volume of the midrapidity region dV/drj = For the parameter set Ng = 0, Nb = 25, 

S/Nb = 200, we choose tm = 5 fm/c which leads to Rt = 5.4 fm. A similar value for tm 
would emerge if an equilibrated QGP with temperature T 500 MeV at the initial time 
Tj = 0.1 — 0.2 fm/c was created ||, |2^. 

For collisions at SPS we employ A^^ = 0 initially. The specihc entropy required to fit the 
particle ratios measured in S' -|- Au with our equation of state is S/Nb = 45 [jT^. For central 
Pb + Pb we employ the same value. In S' -|- Au the net baryon number at midrapidity, 
as measured by the NA35 collaboration 


25| , is dNB/drj ^ 16. 


For Pb + Pb we assume 
dNB/drj = 80 which results in 61 net nucleons at midrapidity (cf. section H). If we 
simply scale by 1/2 (at high energies the initial isospin asymmetry is “deposited” in the 
pions and the nucleons become isospin symmetric), the net proton multiplicity is close to 
the preliminary NA49 data [^. For the (proper) time where the evolution of the mixed 
phase starts we assume tm = 1-5 fm/c (S' -f Au) and tm = 3 fm/c {Pb + Pb), which leads 
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to Rt = 3.7 fm and Rt = 5.9 fm, respectively. 


4 Results for RHIC 


Figure shows the time evolution of A, fig, and fs,QGP = Evaporation 

leads to a faster conversion of QGP into hadronic matter and reduces the lifetime of the 
mixed phase from th — tm = 19 fm/c (only longitudinal boostinvariant expansion) to 
th — tm = 11-6 fm/c. Thus, even if the transverse expansion of the mixed phase is extremely 
slow 1^, the lifetime is considerably shortened by the permanent decoupling of particles close 
to the surface. 

In contrast to the situation for lower vacuum pressure B and initial specihc entropy 
H, 1^, 1^, for S/Nb > 40 and B = 380 MeV/fm^ the temperature in the mixed 
phase (T 160 MeV) is essentially constant and does not increase signihcantly (although 
^QGP^jyQGP increase). The variations of the chemical potentials fig and Pg in time are 
much bigger. Also, strangeness separation between QGP and hadrons occurs {fs,QGP > 0) 
even without evaporation {Ng = 0 in the total system in this case). This is a consequence 
of the fact that from Ng = 0 and Nb > 0 also > 0 follows if hadrons are present (i.e., 
A < 1) P]. For the initial conditions as in figure |l] the evaporation of strangeness (which is 
mainly due to kaons) leads to a net strangeness at the final break-up (i.e. at the time where 
A = 0) of Ng/NB = 11.5%. 

Due to the fact that (for these initial conditions) K, p, A are more abundant in the hadron 
gas than the corresponding antiparticles, they are evaporated earlier and their average freeze- 

As shown in ref. |E3, for 


out times are smaller than those of the antiparticles [pT 
low vacuum pressure B^^^ = 160 MeV (slow hadronization) this phenomenon affects the 
properties of the K — K correlation function. In our model, however, as a consequence of 
the boostinvariant expansion most particles are emitted in the breakup of the purely hadronic 
phase, washing out this signal completely. We obtain equal average freeze-out times for K 
and K\ {Tfo)KT< ~ fm/c (these numbers include the kaons from the final breakup but 
not those from decays of higher-mass resonances). 

The initial conditions in individual events may strongly deviate from our standard set. 
There may well be large fluctuations [^| in A^^ as 'well as in the net baryon number Nb, and 
the entropy per net baryon, S/Nb, if one considers only a small subinterval (e.g. ±0.5 around 
midrapidity) of the total rapidity gap. We assume that even in such events boostinvariance 
around midrapidity is approximately valid. 

In the tables we compile particle and cluster abundancies for our standard initial con¬ 
ditions and a calculation with higher initial net baryon number and an excess of s-quarks. 
The first line gives the sum over evaporation and final breakup. Below, the contribution 
from evaporation only (without final breakup) is quoted. This latter part of antibaryons 
and clusters survives even if the space-time volume of the chemically equilibrated purely 
hadronic phase would be significant (which would lead to cooling below Tc). However, as 
already mentioned in the introduction, this is probably not the case [0, 0, 0- The pion 
multiplicity is perhaps somewhat underpredicted, which may be due to the fact that the 
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K 

K 

N 

N 

A 

A 

830 

100 

92.6 

44.4 

25.8 

13.2 

8.67 

267 

31 

27.1 

11.4 

5.97 

3.19 

2.05 

d 

d 

^He 

^He 

AA 



0.104 

0.0355 

1.20 ■ 10-® 

1.43 ■ 10-^ 

0.0051 



0.024 

0.0066 

2.38 ■ 10-^ 

1.79 ■ 10“® 

0.0010 




Table 1: Multiplicities of hadrons and clusters emitted from the expanding midrapidity cylin¬ 
der. Feeding from decays of higher-mass resonances is included. The hrst line shows the total 
multiplicity (evaporated hadrons plus those from the hnal breakup at time r = th) while the 
second line shows only the evaporated particles. Initial conditions: dS/dr] / dNs/dr] = 200, 
dNs/dr] = 25, dNg/dr] = 0, tm = 5 fm/c {Au + Au at RHIC). 


71 

K 

K 

N 

N 

A 

A 

821 

82.6 

108 

51.5 

20.9 

19.7 

5.49 

265 

25.7 

31.7 

13.2 

4.92 

4.68 

1.32 

d 

d 

^He 

*He 

AA 



0.144 

0.0240 

2.30 ■ 10-® 

6.54-10-^ 

0.0102 



0.033 

0.0045 

4.36 ■ 10-^ 

8.47-10-9 

0.0020 




Table 2: Same as in table 1 but for the initial condition dS/dr] / dNs/dr] = 100, dNs/drj = 
dNs/dr] = 50, tm = 5 fm/c {Au + Au at RHIC). 





effective volume correction should be done different for pions . 

For the standard set of initial conditions (table 1), the higher initial specihc entropy 
(as compared to BNL-AGS and CERN-SPS) leads to considerable antibaryon production 
(iV/iV, A/A > 0.5). In S' + Au at SPS one hnds only N/N 0.1. At RHIC, even the ratio 
of deuterons to antideuterons and helium to antihelium is larger than 0.1. 

For the second set of initial conditions (table 2), of course, the initial values of and 
fiq are considerably higher than in Eg. |I[ The higher pg-values enhance the multiplicity of 
deuterons by 40% and that of by a factor of 2. Also, a strong strangeness separation 
between the QGP and the hadrons in the mixed phase occurs, leading to /s^qgp- values up 
to ~ 1.44. For example, in this calculation increases from 37 MeV initially to 48 MeV 
at the time where the transition to purely hadronic matter is complete. These ps-values are 
comparable to those extracted in ref. for S + Au collisions at SPS (cf. £g. ||), while at 
the same time /Xg is considerably lower. This can be clearly observed in the ratio A/A = 3.6, 
which is larger than N/N = 2.5 {S + Au at SPS: N/N k , 9, while A/A Ri 5). 

If experiments were able to trigger on events (or rapidities) with high net baryon and 
strangeness density, the detection of the so-called MEMOs (meta stable multiple strange 
objects) [^, ^ might be possible. As an example, note that for our second set of initial 
conditions the multiplicity of the lightest MEMO, the AA-cluster (for which we assume a 
mass of mAA = 2200 MeV, i.e. a binding energy of 30 MeV), increases by a factor of 2. A 
similar amount of AA-clusters as in Pb + Pb at SPS is emitted (cf. table 3), however, from 
a larger volume. 

Since at RHIG the initial net baryon number is not very far from zero, even events with an 
excess of antiquarks {Nb < 0) may not be too rare. The initial conditions S = 5000, Ng = 0, 
Nb = —25, tm = 5 fm lead to the same particle and cluster yields as in table 1, except 
that particles and antiparticles have to be interchanged. As a result, the multiplicities of 
antibaryon clusters are enhanced by factors of ~ 3 for antideuterons and ~ 8 for antihelium. 


5 Results for SPS 


As shown in ref. |^, the hadron ratios measured in R-l-AM-collisions at SPS can be htted with 
our equation of state assuming isochronous chemical freeze-out (as in [^) at T = 160 MeV, 
fis = 24 MeV, and fiq = 56 MeV (indicated by the diamond in £g. On the other hand, a 
reasonable fit was also obtained assuming formation of a mixed phase in equilibrium, which 
then hadronizes by evaporating hadrons from the surface, until the whole entropy, baryon 
number etc. is converted into free-streaming hadrons. Due to the neglect of hydrodynamic 
expansion, in this latter model the lifetime of the mixed phase, and thus the time interval in 
which the hadrons freeze out [|I3| , is relatively large (8 — 9 fm/c, cf. £g. In contrast 


to the isochronous freeze-out scenario, the average freeze-out times for the various hadron 
species differ considerably. This is due to the fact that /x^ and /ig are strongly time dependent 
and a transition to purely hadronic matter (which freezes out more or less instantaneously) 
does not take place. 

Here, in our present model (including boostinvariant longitudinal expansion) the lifetime 
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71 

K 

K 

N 

N 

A 

A 

570 

75.9 

53.4 

68.1 

7.48 

16.7 

3.04 

130 

18.1 

10.0 

14.3 

1.00 

2.88 

0.463 

d 

d 

^He 

^He 

AA 



0.373 

4.57-10-^ 

2.44 ■ 10-^ 

3.54- 10-y 

0.0130 



0.0873 

4.20 ■ 10“^ 

8.42 ■ 10-® 

1.71 ■ 10-1° 

2.04- 10-3 




Table 3: Same as in table 1 but for the initial condition dS/dr] / dNs/dr] = 45, dNs/dr] = 80, 
dNs/drj = 0, tm = 3 fm/c {Pb + Pb at SPS). 


of the mixed phase is only 4.3 fm/c, the midrapidity region then enters the purely hadronic 
phase and breaks up. As can be seen in £g. this hnal breakup occurs close to the ‘static-fit’ 
point. Most hadrons are emitted at this point, only few are evaporated^. In the other extreme 
case (only longitudinal expansion without evaporation) the lifetime of the mixed phase for 
our equation of state is th — tm = 5.7 fm/c and the breakup at the time th = 7.2 fm/c 
occurs exactly at the ‘static-fit’ point in £g. Our model thus is close to the sudden freeze- 
out scenario (at some specihc (proper) time) and the freeze-out times of the various hadron 
species do not differ significantly. We repeat that for the considered high specihc entropies 
and bag pressures the temperature is essentially time independent, T = 158 ±2 MeV, while 
the chemical potentials fig and fig vary considerably. 

In Pb + Pb collisions, the lifetime of the mixed phase almost doubles (kg. ^ but still is 
2 — 3 fm/c below the one at RHIC. The resulting hadron multiplicities are shown in table 3. 
The pion and net nucleon multiplicities fit well to the preliminary NA49 data . Assuming 


isospin symmetry for the nucleons, p = N/2, we hnd a ratio p/A = 2.5. Due to the higher 
Pq values, d and ^He are more abundant than at RHIC, while the corresponding antinuclei 
are strongly suppressed. Note that we obtain few less deuterons but many more helium 
nuclei than in the UrQMD model [^. One also observes that the fraction of hadrons that 
are evaporated directly from the mixed phase (i.e. without those emitted at the breakup of 
the purely hadronic phase) is higher at RHIC than at SPS. The longer the lifetime of the 
mixed phase is, the more important evaporation gets. At LHC, a large part of (the entropy 
of) the mixed phase will probably be evaporated (i.e. directly converted into free-streaming 
hadrons) before the huid reaches the purely hadronic phase. 


6 Summary 

In this paper we discussed the evolution of the mixed phase presumably created in heavy-ion 
collisions at SPS and RHIC, employing a model that combines (boostinvariant) longitudinal 
hydrodynamical expansion and particle evaporation off the surface of the huid. For rea- 

^ Thereto^ it is obvious that the (midrapidity) particle ratios in our model are very close to those quoted 
in refs. [O O and we do not need to repeat this discussion here. 
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sonable choices of the initial conditions, our model reproduces the hadron ratios measured 
in S' + Au at SPS (as compiled in |^). The interplay of the evaporation process and the 
hydrodynamical expansion (and vice versa) leads to considerably shorter lifetimes of the 
mixed phase (i.e. to faster hadronization) as compared to the scenarios without evaporation 
(only hydrodynamical expansion I, o, n, n ) or without hydrodynamical expansion (only 
evaporation [^, p!3|). 

We have also calculated the multiplicities of various hadrons, clusters, and hyperon- 
clusters in Pb + Pb at SPS and Au + Au at RHIC (our results for S + Au at SPS are 
not presented here since they are similar to those of ref. [jT^). At SPS, jJLq is large (and in 
particular larger than /i^), implying that nucleon to antinucleon and hyperon to antihyperon 
ratios are big (> 5), in contrast to RHIC, where antinucleons and antihyperons are more 
abundant than at SPS by a factor of 3. 

At RHIC, the initial conditions of the central region may fluctuate. This offers the very 
interesting opportunity to produce clusters of antibaryons and hyperons. In events with 
unusually high net baryon density or net strangeness (at midrapidity) the multiplicity of 
d, and MEMOs may be enhanced by up to 100% (as compared to the average). In 
this case > fiq can be reached, leading to an emission of more hyperons per antihyperon 
than of nucleons per antinucleon. On the other hand, negative net baryon number (excess 
of antiquarks) at midrapidity increases antinuclei abundancies by up to a factor of 10. 


Acknowledgements: We thank D. Ardouin and D.H. Rischke for helpful discussions, and 
K.J. Eskola for clearifying comments on entropy production in the preequilibrium stage. 
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Figure 1: Time evolution of the volume fraction occupied by QGP, the chemical potentials 
of s- and u, d-quarks, and the number of net strange quarks per baryon in the QGP. The 
initial conditions are as in table 1 (average Au + Au at RHIG; bottom: only longitudinal 
hydrodynamic expansion without particle evaporation; top; including evaporation). 


12 








^ (MeV) 


Figure 2: Time evolution of /i^ and fig in the evaporation model of ref. (without boost- 
invariant hydrodynamic expansion), compared to our model. Our initial conditions for the 
midrapidity cylinder are dS/drj / dNs/dr] = 45, dNs/dr] = 16, dNs/drj = 0, tm = 1-5 fm/c 
(S' -|- Au at SPS). The parameters of the static fit are also indicated jl^ . 
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Figure 3: Time evolution of Hs and jiq in our model (including evaporation and boostinvariant 
expansion). Initial conditions for the midrapidity cylinder as in table 3 {Pb + Pb at SPS). 
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